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Some phase transitions of cosmological interest may be weakly first-order and cannot be analyzed by 
a simple perturbative expansion around mean field theory. We propose a simple two-scalar model — 
the cubic anisotropy model — as a foil for theoretical techniques to study such transitions, and we 
review its similarities and dissimilarities to the electroweak phase transition in the early universe. 
We present numerical Monte Carlo results for various discontinuities across very weakly first-order 
transitions in this model and, as an example, compare them to e-expansion results. For this purpose, 
we have computed through next-to-next-to-leading order in e. 



Renormalization-group (RG) techniques for studying 
phase transitions, such as the e expansion, have been 
used for over twenty years, and the analysis of second- 
order transitions has produced quite successful quanti- 
tative predictions Weakly first-order transitions, 
in which the transition dynamics is dominated by long 
but finite wavelength fluctuations, were also investigated 
using e expansion techniques However, there has 
been less emphasis on testing quantitative predictions for 
such first-order transitions, in part because few predic- 
tions were easily accessible to experiments of the time [jfj . 

In recent years, renewed interest in quantitative pre- 
dictions for weakly first-order transitions has arisen in 
cosmology, specifically the genesis of matter. Current 
electroweak theory holds the promise of explaining the 
net baryon number of the universe Any scenario 
of baryogenesis must satisfy Sakharov's conditions of (i) 
baryon number violation, (ii) C and CP violation, and 
(Hi) departure from thermal equilibrium. Electroweak 
theory has all of these ingredients: baryon number viola- 
tion comes from a quantum anomaly, C and CP violation 
is built in, and a first-order electroweak phase transition 
can provide the departure from equilibrium. The elec- 
troweak transition separates the low-temperature world, 
where the SU(2)xU(I) gauge symmetry is hidden and 
only the U(l) symmetry of electromagnetism is mani- 
fest, from the high-temperature domain where the full 
SU(2)xU(l) symmetry is manifest. To test the viabil- 
ity of electroweak baryogenesis, detailed understanding 
of many properties of this phase transition is required. 

Fig. |l| sketches what can be firmly established about 
the phase diagram of electroweak theory in the mini- 
mal Standard Model (the standard model with a single 
SU(2) doublet Higgs field) without resorting to numerical 
simulations. A more complicated Higgs sector is prob- 
ably needed for realistic electroweak baryogenesis, but 
the minimal case is easier to discuss and provides a use- 
ful testing ground. The one important parameter of the 
minimal model not well determined by experiment is the 
Higgs boson mass mi,, or equivalently, the Higgs quartic 
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FIG. 1. Phase diagram of electroweak theory, or the cubic 
anisotropy model (for u < 0,v > 0), based purely on analytic 
arguments. The left and right axes show the relevant ratio R 
of couplings for electroweak theory, or for the cubic anisotropy 
model, respectively. 



self-coupling A. For a sufficiently light Higgs mass (or 
small A) the phase transition can be successfully treated 
by expanding around mean field theory and is known 
to be fluctuation-induced first order. In other words, it 
appears to be second-order in mean field theory, but per- 
turbative fluctuations generate a finite correlation length 
and render the transition first order. For large Higgs 
mass (or large A), the problem becomes difficult to an- 
alyze. Perturbative techniques are only reliable when 
mi,<Mw, where M-yv ~ 80 GeV is the W-boson mass. 
In condensed matter applications, this is commonly re- 
ferred to as the Ginsburg criteria for the success of mean 
field theory. An equivalent way of stating this criteria is 
that the ratio of couplings R = X/g^,, where g^ is the 
electroweak gauge coupling, needs to be small. If one 
(unrealistically) sends g w (and hence M\v) to zero, so as 
to probe the opposite R — > oo extreme, there is one piece 
of solid information. The Higgs sector decouples from 
the gauge interactions in this limit, and becomes a pure 
0(4) symmetric scalar theory, which is known to have a 
second-order transition. 

Viable electroweak baryogenesis generally needs as 
strong a first-order phase transition as possible, and so 
prefers as small a Higgs mass as possible. Experiment, 
however, puts lower bounds on Higgs masses (60 GeV 
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in the minimal model) and forces one into the regime 
where a perturbative analysis of the transition is at best 
marginal. Consequently, there is a need to understand 
the full phase diagram for electroweak theory and to 
develop techniques for reliable computations even when 
?7ih ^ Mw 

There arc several things known about the transition in 
various unphysical limits. Near but just below 4 spatial 
dimensions (d=4 — e), the transition remains first order 
for large Higgs mass ||,[L0]], as we review below. In the 
large N (scalar) limit (or in 2 + e dimensions), there is 
a value of rrih at which there is a tricritical point, and 
the transition becomes second-order for larger mi, jl0|,[nj . 
In contrast, recent lattice simulations suggest that the 
first-order line actually ends in a critical point at mh ~ 
above which there is no phase transition whatsoever 
fiJI . There has been a recent attempt to understand this 
behavior using RG methods in three dimensions ||l3|| . 

In all the cases outlined above, one is left with the prob- 
lem of how to analyze the phase transition for those Higgs 
masses where it is only weakly first-order and where per- 
turbation theory breaks down. For any techniques that 
purport to address this problem, it would be useful to 
have an extremely simple statistical system — on the level 
of simplicity of the Ising model — that could serve as a 
canonical initial testing ground. The Ising model pro- 
vides a canonical example of a second-order transition 
and is in the same universality class as the continuum 
field theory of a single scalar field. To study weak first- 
order transitions, one may add a second spin field to make 
a model known as the Ashkin- Teller model, or equiva- 
lently add a second scalar field to make a model known 
as the cubic anisotropy model. 

The purpose of the present work is to review these 
models; review the similarities and dissimilarities with 
the electroweak phase transition (see also present 
the results of numerical simulations of first-order transi- 
tions; and, as an example of testing a theoretical method, 
compare the numerical results to predictions of the e ex- 
pansion. The e expansion is the generalization of a theory 
from three to 4— e spatial dimensions, where the dynam- 
ics of long wavelength fluctuations is tractable using RG- 
improved perturbation theory, followed by extrapolation 
of results back to e—1. We have studied weakly first- 
order transitions of the cubic anisotropy model through 
next-to-next-to-leading order in e. (Leading order re- 
sults were obtained twenty years ago in ref. j|.) We 
propose the cubic anisotropy model as a benchmark for 
other investigators who advocate particular techniques 
for studying weakly first-order transitions. The details 
of our e expansion calculations and numerical studies are 
given elsewhere 

In its general form, the cubic anisotropy model is an 
O(n) symmetric scalar model to which an additional in- 
teraction is added which breaks O(n) down to hyper- 
cubic symmetry ||: 
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We shall focus on the simplest case, n=2. This model is 
analogous to electroweak theory when u < and v > 0. 
Fig. |l] illustrates the phase structure for weak couplings 
(|u|, v <C 1). There is again a suitable ratio of couplings, 



R = (v — n\u\ 
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which defines a Ginsburg criteria. When the ratio drops 
below zero, the tree-level potential in ([!]) is unstable, just 
as A < makes the Higgs potential unstable in elec- 
troweak theory. When the ratio is small, perturbation 
theory turns out to be well behaved. Then the phase 
transition is first order and, like electroweak theory, is 
fluctuation driven. When the ratio is moderate or large, 
perturbation theory breaks down. When it is infinite, 
corresponding to u=0, the theory reduces to n decoupled 
copies of a theory of a single scalar field. Each such copy 
is in the same universality class as the Ising model, and 
has a second-order transition. 

To further appreciate the similarities between the cu- 
bic anisotropy model and the electroweak phase transi- 
tion, we turn to the e expansion. When e is small, both 
theories can be solved using RG-improved perturbation 
theory. Fig. ||a depicts the RG flow of the couplings g w 
and A of electroweak theory in 4— e spatial dimensions 
as one examines increasingly large distance scales. Not 
shown is the third relevant parameter, the running Higgs 
mass, which always grows (relative to the renormaliza- 
tion point) as one scales to larger distances. For g w ^0, 
the couplings do not flow to any fixed point. They do, 
however, flow into the region A <C g^, where the Gins- 
burg criteria is satisfied. After using the renormalization 
group to flow into this region, one can use standard per- 
turbative techniques to make quantitative computations 
of thermal properties and conclude that the transition 
is first order. This is the argument that, for sufficiently 
small e, the line of first order transitions shown in Fig. [l] 
continues all the way to R — oo. 

In condensed matter physics, the underlying theory at 
short distances is typically strongly interacting. In con- 
trast, the Higgs model in particle physics is generally of 
interest as a fundamentally weakly coupled (X,g w *C 1) 
theory at short distances |17|j . Hence, one can focus 
attention on short-distance theories defined by varying 
A/g^ in the neighborhood of the Gaussian fixed point 
at A=g w =0, as illustrated by the arc A near the origin 
in Fig. 0a. By renormalization group equivalence, this 
covers all theories in the lightly shaded region of fig. ||a. 
In particular, it covers the case of the real world (g w 
fixed but A unknown, designated by line J3), provided 
the Higgs mass is not so large that the scalar sector is 
fundamentally non-perturbative (i.e., mi, < 1 TeV). 

Analogous flows for the cubic anisotropy model are 
show in fig. 0b. For u < and v > 0, the flows are 
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FIG. 2. Renormalization group flows, for small e, of (a) 
electroweak theory, and (b) the cubic anisotropy model. 

qualitatively similar to the electroweak case. For u = 0, 
the theory flows to the Ising fixed point. For u small 
and negative, it instead flows to a region where the Gins- 
burg criteria is met and the transition is known to be first 
order — an observation first exploited by Rudnick Q . The 
Ising fixed point is a tricritical point separating first-order 
(u < 0) and second-order (u > 0) transitions. 

An important distinction between the cubic anisotropy 
model and electroweak theory is that the cubic 
anisotropy model must have a phase transition for any 
choice of the Ginsburg ratio R because a global symme- 
try is spontaneously broken. No such argument holds for 
electroweak theory due to Elitzur's theorem fl§|| . 

In analogy with electroweak theory, one could attempt 
to analyze the phase transitions in the cubic anisotropy 
model for weak couplings u and v, and for arbitrary 
choices of the Ginsburg ratio R. However, we know that 
the transition gets weaker, and the analysis will become 
harder, as R increases. So an instructive goal is to ana- 
lyze the extreme case R — > oo. The limiting renormaliza- 
tion group trajectory is the one marked C in Fig. ||b 
that runs from the Gaussian fixed point to the Ising 
fixed point, and then away from the Ising fixed point 



to the perturbative regime. This extreme case has a sub- 
stantial advantage for numerical simulations: the same 
long-distance theory in the perturbative regime can be 
reached starting from a strongly coupled short-distance 
theory described by the trajectory D which comes in from 
(u=0, v=oo) to the Ising fixed point and then flows away 
to the perturbative regime along C. The advantage of 
starting with this strongly coupled theory is that it can be 
replaced by a discrete spin model, which is easier to sim- 
ulate numerically than theories with continuous-valued 
fields. 

In the (lattice regularized) theory of a single scalar field 
with 

S = J d d x [\{dcpf + A(</> 2 - i; 2 ) 2 ] , (3) 

the field 4> is constrained to only take on the two values 
±u when A — > oo, and hence the theory reduces to a 
simple Ising model with a discrete two component spin. 
Similarly, the strong coupling limit of the n = 2 cubic 
anisotropy model is a model of two coupled Ising spins s 
and t known as the Ashkin- Teller model [l9|j20[| : 

S = -py^^SjSj + Utj + xSiSjtit-j) , (4) 

m 

where the sum is over all neighboring pairs of sites (i, j); 
the spins s,t = ±1; and j3 and x are parameters that rep- 
resent the inverse temperature and the deviation from 
the limit of decoupled Ising models, respectively. The 
R — > oo limit in our discussion of the cubic anisotropy 
model corresponds to the x^Q + limit of the Ashkin- 
Teller model. 

The x-^0 + limit is crucial if one wants the Ashkin- 
Teller model to correspond to the continuum cubic 
anisotropy model. That is because generic first-order 
transitions have finite correlation lengths of 0(1) lattice 
spacings at the transition. Everything about such transi- 
tions depends on the details of the short-distance physics: 
the type of lattice used, the exact details of the cou- 
plings, etc. In the x^0 + limit, however, the correlation 
length at the transition becomes arbitrarily large (since 
for x — it is infinite), and so in this limit the details of 
the short-distance physics decouple. Dimensionful quan- 
tities like the correlation length £ aren't directly useful 
quantities to predict in this limit, since they diverge, but 
ratios such as 

ton ft = hm p^iZ i (5) 

are non-trivial and universal. Here £ + is the correlation 
length infinitesimally above the critical temperature, and 
£_ is that infinitesimally below. Our project has been 
to compute such ratios in the e expansion for the cubic 
anisotropy model and to compare to numerical simula- 
tions of the corresponding Ashkin- Teller model. 
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The details of our e expansion calculations are given in 
ref. fuJl. We find 



— = 2 [l + e (-0.1063 lne + 0.4494) 
X- 

+ e 2 (-0.00731n 2 e + 0.1647 lne - 0.2859) + 0(e 3 



(6) 



e 

c 
c_ 



2 [l + e (-0.1063 lne + 0.4139) + 0(e 2 



(7) 



7T- = m I e + e 2 (-0.21251ne + 1.3993) + 0(e 3 )]. (8) 



C is the specific heat, x = d{4>i) / dh\h=a is the scalar 
susceptibility, defined by adding an interaction h<j>\ to 
the model, and to leading order is the same as £ 2 . Our 
leading-order results for the \ an d C ratios differ by fac- 
tors of 4 from those originally reported by Rudnick M . 

Table | summarizes the leading-order and next-to- 
lcading order results for our three ratios. In general, 
e expansions are asymptotic expansions and the coeffi- 
cients begin to grow after just a few terms. For critical 
exponents, one can often get reasonably good values by 
simply adding terms at e=l until the contributions start 
to increase. Our results in the cubic anisotropy model 
are not as well-behaved as the best Ising model series, 
but are comparable to the e expansion for the exponent 
rj. The expansion for \+ /x— leads one to expect that the 
actual result should lie somewhere between the LO and 
NLO results. One might have a similar expectation for 
C-|_/C_ and even though the expansion for the 

former is clearly poorly behaved. 

Very accurate predictions for critical exponents in the 
Ising model have been obtained using resummation tech- 
niques which combine knowledge of the asymptotic large 
order behavior of the series with explicit low-order coeffi- 
cients Our series expansions, however, contain loga- 
rithms of e, and we do not currently know how to derive 
the large-order behavior of these series. As explained in 
ref. |l5j, these logarithms arise because physics at the 
phase transition involves two scales whose ratio is O(e). 

Our numerical results are also presented in Table B, and 
details may be found in ref. [jl6| . One sees that C + / C_ is 
indeed bracketed by the LO and NLO e-expansion results 
and that £+/£- works reasonably well, x+lx-i however, 
is a disappointment for the e expansion, since the numer- 
ical result is not bracketed by LO and NLO, and the the 
NNLO correction only makes the agreement worse. 

Our results highlight the utility of the cubic anisotropy 
and Ashkin- Teller models as foils for theoretical methods 
for making quantitative predictions of weakly first-order 
transitions. It would be interesting to see the comparison 
of our numerical results against methods other than the e 
expansion, and it would be useful to have more accurate 
Monte Carlo determinations of the universal ratios. 

This work was supported by the U.S. Department 
of Energy grants DE-FG06-91ER40614 and DE-FG03- 
96ER40956. 



universal ratio 




e expansion 




Monte Carlo 




LO 


NLO 


NNLO 




C+/C- 


0.053 


0.127 




0.069(8) 


t+li- 


1.4 


1.7 




1.7(2) 


x+/x- 


2.0 


2.9 


2.3 


4.1(5) (a) 



TABLE I. Leading (LO), next-to-leading (NLO), and 
next-to-next-to-leading order (NNLO) results in the e expan- 
sion compared to numerical Monte Carlo results. The error 
estimate on the Monte Carlo measurement of x+/x- should 
be taken with a grain of salt (see ref. lid] for details). 
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